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^ The common Muon and Proton Apparatus for Structure and Spectroscopy (COMPASS) is a 

multi-purpose fixed-target experiment at the CERN Super Proton Synchrotron (SPS) aimed at 
I studying the structure and spectrum of hadrons. 

In the naive Constituent Quark Model (CQM) mesons are bound states of quarks and antiquarks. 

QCD, however, predict the existence of hadrons beyond the CQM with exotic properties inter- 

^ preted as excited glue (hybrids) or even pure gluonic bound states (glueballs). One main goal of 

1) COMPASS is to search for these states. Particularly interesting are so called spin-exotic mesons 

pi I which have J^^ quantum numbers forbidden for ordinary qq states. 

^ Its large acceptance, high resolution, and high-rate capability make the COMPASS experiment an 

excellent device to study the spectrum of light-quark mesons in diffractive and central production 
reactions up to masses of about 2.5 GeV/c'. COMPASS is able to measure final states with 
charged as well as neutral particles, so that resonances can be studied in different reactions and 
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decay channels. 

During 2008 and 2009 COMPASS acquired large data samples using negative and positive sec- 
ondary hadron beams on £H2, Ni, and Pb targets. The presented overview of the first results from 
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this data set focuses in particular on the search for spin-exotic mesons in diffractively produced 
n^n^n^ , rjn, rj'n, and n^n'^n^n'^n^ final states and the analysis of central-production of 
TT+TT" pairs in order to study glueball candidates in the scalar sector. 
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1. Experimental Setup 

COMPASS is a two-stage high-resolution spectrometer that covers a wide range of scattering 
angles and particle momenta [1]. Both stages of the spectrometer are equipped with hadronic 
and electromagnetic calorimeters so that final states with charged as well as neutral particles can 
be reconstructed. A Ring-Imaging Cherenkov Detector (RICH) in the first stage can be used for 
particle identification. It is able to separate kaons from pions up to momenta of about 50 GeV/ c. 
The target is surrounded by a Recoil Proton Detector (RPD) that measures the time of flight of recoil 
protons using two scintillator barrels. COMPASS uses the M2 beam line of the SPS which can 
deliver secondary hadron beams with a momentum of up to 280 GeV/c and a maximum intensity 
of 5 • 10^ s^^ The negative hadron beam that was used for the analyses presented here has a 
momentum of 190 GeV/c and consists of 97 % n^,2 % K , and 1 % ^ at the COMPASS target. 
Two ChErenkov Differential counters with Achromatic Ring focus (CEDAR) upstream of the target 
are used to identify the incoming beam particles. 



2. Search for Spin-Exotic Mesons in Diffractive Dissociation 

Diffractive dissociation reactions are known to exhibit a rich spectrum of produced intermedi- 
ate states. In the past several candidates for spin-exotic mesons have been reported in pion-induced 
diffraction [2]. In diffractive events the beam hadron is excited to some intermediate state X via 
^-channel Reggeon exchange with the target. At 190 GeV/c beam momentum Pomeron exchange 
is dominant. The X decays into a n-body final state which is detected by the spectrometer. The 
process beam + target — )• X + recoil, where X — )• /ii . . . /j„, is characterized by two kinematic vari- 
ables: the square of the total center-of-mass energy, s, and the squared four-momentum transfer to 
the target, t = (pbeam — Px)~- It is customary to use the variable ?' = |f| — 1?^;^^ instead of t, where 
\t\^^ is the minimum value of \t\ for a given invariant mass of X. 

2.1 K'Tf^Tt' Final States from K~ Diffraction 

A partial- wave analysis (PWA) of n^Ti^n^ final states produced in diffraction on a Pb 
target in the squared four-momentum transfer range of 0.1 <t' < 1.0 (GeV/c)^ and extracted from 
data taken during the pilot run in 2004 showed significant intensity in the spin-exotic J^'^ = 1^+ 
partial wave in the p"?!^ decay channel [3]. A Breit-Wigner fit yielded resonance parameters 
consistent with the disputed tti (1600) claimed in this channel by other experiments [4]. 

In 2008 COMPASS has acquired a large data set of the diffractive dissociation reaction p — )• 
71^71+ 71^ Psiow The trigger included a beam definition and the RPD, which ensured that the target 
proton stayed intact and also introduced a lower bound for t' of about 0.1 (GeV/c)^. Diffractive 
events were enriched by an exclusivity cut around the nominal beam energy. After all cuts the 
analyzed n^n^n^ sample contains about 5 • 10^ events in the range 0.1 < < 1.0 (GeV/c)^. 

The PWA approach employs the isobar model [5] in which the decay X — )• n^n^n^ is de- 
composed into a chain of successive two-body decays. The X with quantum numbers Z^*- and 
spin projection is assumed to decay into a ;r+7r^ resonance, the so-called isobar, and a bache- 
lor pion. The isobar has spin S and a relative orbital angular momentum L with respect to ^^[^(.heior- 
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Figure 1: Left and center columns: n^n^n^ invariant mass spectrum (top left). Intensities of the three 
major waves: 1++0+ [p7i\S wave with the ai(1260) (top center), 2++ 1+ [p7r]£> wave with the a2(1320) 
(bottom center), and 2^+0+[/27r]5 wave with the 7^2(1670) (bottom left). Right column: Examples for 
small waves (note the different y scale): 0^+0+ [/o(980);r]5 wave with the ;r(1800) (top right) and 2++ 2+ 
[pn\D with the 02(1320) (bottom right). 



A partial wave is thus defined by 7^''M'^[isobar]L, where e = ±1 is the reflectivity which corre- 
sponds to the naturality of the exchanged particle in the production process [6]. 

The spin-density matrix for the chosen set of 52 partial waves plus an incoherent isotropic 
background wave is determined by unbinned extended maximum likelihood fits performed in 
20 MeV/c^ wide bins of the three-pion invariant mass mx- In these fits no assumption is made 
on the produced resonances X other than that their production strength is constant within a mx 
bin. The PWA model includes five n+n- isobars [7]: (7r7r)5.wave, P(V70), /o(980), /2(1270), 
and p3(1690). They were described using relativistic Breit-Wigner line shape functions including 
Blatt-Weisskopf barrier penetration factors [8]. For the 71+ TT^ 5-wave we use the parametriza- 
tion from [9] with the /o(980) subtracted from the elastic nn amplitude and added as a separate 
Breit-Wigner resonance. Mostly natural-parity waves are needed to describe the data. A rank-two 
spin-density matrix was used in order to account for spin-flip and spin-non-flip amplitudes at the 
target vertex. 

The intensities of the three dominant waves in the n^n^n^ final state, 1++0+ [p7l]S, 2++ 1+ 
[p7l]D, and 2^+ 0+ [f27l]S, are shown in Fig. 1 left and center. They contain resonant structures that 
correspond to the well-known ai (1260), (32(1320), and 712(1670), respectively [7]. That the applied 
analysis technique is able to cleanly separate even small waves with intensities at the percent level 
of that of the dominant waves is illustrated in Fig. 1 right. It shows two exemplary waves, 2++ 2+ 
[pK]D and 0^+0+ [/o(980)7r]S, which exhibit clear peaks of the a2(1320) (this time with M = 2) 
and the 7r(1800), respectively. 

Figure 2 left shows the intensity of the spin-exotic 1^+1+ [p7i]P wave. The bump around 
1.2 GeV/c^ does not seem to be of resonant nature. It is unstable with respect to changes in the 
PWA model which hints that it is rather an artifact of the analysis method. On the other hand the 
peak structure around 1.6 GeV/c^ as well as the corresponding rising phase motion with respect to 
the tail of the a\ (1260) in the 1++ 0+ [p7i]S wave (cf. Fig. 2 center) are stable against modifications 
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Figure 2: Le/f; Intensity of the spin-exotic 1 +1+ [pTtjD wave. Center: Phase of the 1 + wave relative to 
the 1++ 0+ \pK\S wave. Right: Relative phase with respect to the 2^+0+ [/2?r]5 wave. 



of the fit model. As Fig. 2 right shows, the structure is phase locked with the 7r2(1670) in the 2^+ 
0+ [/27r]5' wave. This is consistent with the results obtained from a PWA of the pilot-run data taken 
with a Pb target [3]. The interpretation of the 1^+ wave in terms of resonances, however, is still 
unclear. There seem to be significant contributions from non-resonant Deck-like processes that 
need to be included into the fit model. 

The partial waves also exhibit different t' behavior depending on their spin projection M and 
their resonance contents. This is problematic, since the data cover a rather wide t' range from 0. 1 to 
1.0 {G&W/cY while at the same time the PWA model assumes full coherence of all partial waves. 
This issue will be addressed by performing a two-dimensional PWA in bins of nix and t' . 

2.2 ic~r\ and K~r\' Final States from K~ Diffraction 

Previous experiments claimed spin-exotic J^^ = 1^+ resonances also in the n^r] [10, 11] 
and n^f]' [10, 12] final states. However, the resonant nature of the observed signals is still 
controversial [13]. COMPASS performed a PWA of both final states in the diffractive reaction 
p —?■ ?r"T]^''psiow TT^Tr+TT" 77Psiow in the kinematic range 0.1 < ?' < 1.0 (GeV/c)^ [14]. 
The T] is reconstructed via its decay to 7r+7r^7r" where the 71° goes to a photon pair. In order to 
reconstruct the Tj', T] reconstructed in the 77 channel are combined with a n^Ti^ pair. The final 
selected data sample contains about 100000 n^rj and 35 000 tt^t]' events. 

The performed PWA follows previous analyses and includes S, P, and D waves with M < I 
and both natural and unnatural parity exchange. In addition a 4++ 1+ G wave and an incoherent 
isotropic background wave were included. For both channels the bulk of the intensity is described 
by the 1^+, 2++, and 4++ waves with M = I and natural parity exchange. 

The left column of Fig. 3 shows the 2++ 1+ (D+) intensity for the two final states. Although 
the two distributions look rather different they are related to each other by the fact that r] and 17' 
are mixtures of the basis states T]„ = ( 1 / ^/2) {uu + dd) and r]s = ss. The partial-wave amplitudes 
Tj of a spin-/ resonance with negligible ss content decaying into n^r] and n^r]' should be related 
to each other by the T]-T]' pseudoscalar mixing angle in the flavor basis, phase space, and barrier 
penetration factors [15]. With q being the two-body breakup momentum the amplitude ratio can be 
written as 

ffn' r / "1 /+1 /2 



tan0 



(2.1) 



Here simplified baiTier factors of the form q-' were used. 

The bottom left plot in Fig. 3 shows in red the intensity of the rj channel scaled by the 
phase space factor from Eq. (2.1) overlaid on top of the corresponding intensity in the n^rj' final 
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Figure 3: Left column: Intensity of the 2++ 1+ (£>+) wave in the n^rj (top) and n^rj' (bottom, black) 
channel, respectively. The red dots in the bottom plot show the rj D+ intensity scaled by phase space. 
Center and right column: Intensities and relative phases for tt^tj' (black) and n^i] scaled by phase space 
(red): 4++1+ (G+) intensity (top center), spin-exotic 1^+1+ {P+) intensity (top right), relative phases 
G+ — D+ (bottom center) and D+ — P+ (bottom right). 



state. The phase space scahng leads to a remarkable similarity of the n^r\ and tt^T]' intensities. 
This behavior is expected for light-quark resonances like the a2(1320), but it is rather astonishing 
that the scaling also holds in the high-mass region where non-resonant contributions are expected 
to dominate. As the middle row of Fig. 3 shows, the same behavior is observed for the 4++ 1+ 
(G+) partial waves with the 04(2040). Also the relative phase between the G+ and the waves 
is very similar in both channels which means that the physical composition of the two waves is 
almost identical. 

The picture is much different for the spin-exotic 1^+1+ (P-^_) wave which is the dominant wave 
in the n^r\' channel with a prominent broad structure around 1.6 GeV/c^ whereas it is strongly 
suppressed in the 7i^r\ channel (cf. Fig. 3 right). This observation is consistent with previous 
experiments and with the suspected non-^^ character of this wave [16]. The P+ wave exhibits a 
slow phase motion with respect to the wave in the 1.6 GeV/c^ mass region which, however, 
evolves differently in the two final states. This is a hint that the resonant content is different. The 
resonance interpretation of the /+ wave requires a better understanding of the high-mass region of 
the D+ and G+ waves against which the phase motion of the Pj^ is measured. In this region non- 
resonant contributions from double-Regge processes are expected to play an important role which 
should be included in the fit model [17]. 

2.3 1t~1t^1t~1t^1t~ Final States from lt~ Diffraction 

The diffractive reaction n^Vh — 71^71+71^71+71^ Pb is very interesting, because the 57r final 
state gives access to the mass region around and above 2 GeV/c^ which is also called the "light- 
meson frontier". Only little is known about resonances in this region. There are many missing 
states as well as states that need confirmation. 

A PWA of 200 000 exclusive events was performed in the kinematic region < 5 • 10^^ (GeV/ c)^. 
Since a Pb target was used the recoil could not be detected. The partial-wave decomposition of 
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TC'^TC Isobars 



Tt n'^n^ Isobars n^n n^n Isobars 



(?r7r)5-wave,p(770) flj (1260), fl2(1320) /2(1270),/i (1285), /o(1370, 1500), p'(1450, 1700) 
Table 1: Isobars used in the wave set of the n^n'^n^n^n^ PWA [18]. 



diffractively produced 7r^7r+7r^7r+7r^ final states is cliallenging. Due to tlie lack of dominant 
structures in tlie 5 71 invariant mass distribution and that of its subsystems, it is difficult to con- 
straint the wave set which is potentially large, because of the many allowed isobars. In order 
to explore the model space more systematically an evolution- 
ary algorithm was developed that uses a Bayesian goodness- 
of-fit criterion that takes into account the model complexity 
to iteratively find the best set of waves [18]. A pool of 284 
waves was offered to the algorithm. The best model found 
consists of 3 1 partial wave plus an incoherent isotropic back- 
ground wave. Table 1 lists the used isobars. An advantage 
of this automatized method is that it not only finds a wave 
set but also gives an estimate of the systematic uncertainty 
introduced by the particular choice of the wave set. 

The mass dependence of a submatrix of the spin-density 
matrix consisting of 10 out of the 31 waves was parametrized 
using a simple model that consists of 6 resonances that were 
described by constant-width relativistic Breit-Wigners plus 
background terms. Although mixing and coupled-channel 
effects were neglected the model is able describe the data 
surprisingly well as is illustrated in Fig. 4. Table 2 sum- 
marizes the extracted resonance parameters. Being the first 
analysis of diffractively produced , the re- 

sults of this PWA should be interpreted with care. The anal- 
ysis is still based on a number of assumptions. In particular 
the lack of knowledge on the properties of the four-pion iso- 
bars is difficult to quantify. 



Resonance 




Fit result 
MeV/c^ 
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1400 (at Umit) 
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500 (fixed) 
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1781 ±5+^ 
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443 ± 14+^^ 
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402 ± 17+^f 


7^2(1670) 
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1719.0 (fixed) 
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251.4 (fixed) 


712(1880) 
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1854 ±6+^ 
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259±13;^i 


712(2100) 


M 


2133 ±12^^^^ 




r 


448±22;5[J 



Table 2: Summary of extracted 57r 
resonance parameters [18]. The first 
uncertainty is statistical, the sec- 
ond represents the systematic er- 
ror including uncertainties from the 
choice of the final wave set. 



3. Search for Scalar Glueball Candidates in Central Production 



Pomeron-Pomeron fusion processes are believed to 
provide a glue-rich environment which should lead to an 
enhanced cross section for the production of glueballs. 
Although glueball candidates are discussed in the litera- 
ture, their existence could not yet be confirmed experi- 
mentally. Continuing the efforts of CERN OMEGA spec- 
trometer in the late 1990s [19] COMPASS studies the 
central-production reaction pp — )• Piast^"^^^ Psiow using a 
190 GeV/c positive secondary hadron beam which consists 
of 75 % p, 24 % 71+ , and 1 % K+ at the ^2 target [21]. 




Figure 5: 7r+7r invariant mass dis- 
tribution. 
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Figure 4: Fit of the m^jc dependence of a spin-density submatrix of 10 waves [18]. The diagonal elements 
show the wave intensities, the off-diagonal ones the interference terms (top right triangle: real, bottom left 
triangle: imaginary part). The light blue error bars represent the uncertainty from the choice of the wave set. 



Centrally produced tt+tt" are separated from beam- 
diffraction events by cutting on the invariant mass m(/7fast7r^) > 1.5 GeV/c^. Elastic scattering 
at the target vertex is ensured by the RPD trigger. This selects tt+tt^ pairs within \xf\ < 0.25. At 
COMPASS energies in addition to Pomeron-Pomeron fusion also Reggeon-Pomeron and Reggeon- 
Reggeon processes are expected to contribute. This can be seen in the 71+ 71^ invariant mass spec- 
trum shown in Fig. 5, where in addition to structures from /o(980) and /2(1270) a clear p(770) 
peak is visible. 

The employed PWA method is similar to the one used by WA102 [20]. The n'^n" pair is 
assumed to be produced in the collision of two exchanged objects emitted by the scattered protons. 
The one emitted by the beam proton caiTies a squared four-momentum transfer ti, the one from 
the target proton t2, respectively. We make the strong assumption that exchange particle ti only 
transmits helicity A = and that we can treat it like an external particle. 

A complication arises from the fact that the two-pseudoscalar final state suffers from mathe- 
matical ambiguities, which means that for a given wave set different partial-wave decompositions 
result in exactly the same angular distribution. The PWA was performed using a wave set con- 
sisting of S, P, and D waves which leads to eight ambiguous solutions. Additional constraints are 
needed to select the physical solution. For six solutions, most of the intensity accumulates in a 
single wave which is clearly unphysical, since at least one resonance is known to be present in S, 
P and D waves (cf. Fig. 5). The choice among the two remaining solutions is not evident, the 
intensities and phases are very similar. The remaining solutions are compatible with the physical 
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constraints. A clear /2(1270) peak can be seen in the 2++0~ (Z)q ) intensity. It is accompanied 
by a phase motion with respect to the O^+O" (5q ) wave. The P waves show p(770) peaks. The 
most interesting wave is, of course, the Sq . Its interpretation is, however, challenging and work in 
progress. 

4. Summary 

COMPASS is a unique apparatus to study light-quark hadron spectroscopy in diffractive and 
central production reactions. Large data sets were collected for various final states with charged as 
well as neutral particles using different target materials. The main focus of the first analyses lies 
on the search for spin-exotic mesons. J^^ =1"+ signals are observed in diffractively produced 37i 
and KT] ^'^ final states. However, their interpretation in terms of resonances is complicated by strong 
non-resonant contributions that need to be taken into account in the analysis. In the future, the 
search for spin-exotic signals will be extended to channels like n^con^^ and {KK)^7i^ . Glueball 
candidates are studied in centrally produced n^n^ pairs. First PWA show promising results and 
this analysis is beeing extended to K^K~ , k^k^, and K'^^l final states. 
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